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Abstract 

We study the propagation of waves in a quasi ID homogeneous host medium filled with various res- 
onators. We first prove that a far field coupling between the elements explains its dispersive nature. This 
coupling is interpreted as a Fano interference between the incoming wave and the waves re-radiated by 
each resonator, which experience a phase shift at resonance. We propose a simple formalism that gives 
the complete dispersion relation of the medium in terms of the far field response of a single resonator. We 
prove that our approach applies to and unifies various domains such as metamaterials, hybridization band 
gap materials and designer's plasmons. Finally we show that those media, spatially random or organized 
on a scale larger than the wavelength, also present very interesting properties, which broadens the range of 
man made exotic materials. 
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Propagation media containing resonant inclusions have been studied for over a century in 
acoustics, electromagnetism or solid state physics. There exist some in nature, such as dielectrics, 
which contain enormous amounts of atoms. To calculate those materials permittivities one con- 
siders that each atom "sees" the same electromagnetic field and calculates the average field that 
takes into account an incoming wave as well as the overall response of the ensemble of atoms 
ifTI . This macroscopic view assumes that there are no variations of the electromagnetic field at the 
scale of the inter-atomic distance. Later on, the concept of polariton has been proposed to explain 
the propagation of electromagnetic waves in dielectrics [|2l[3]|. The polariton approach of the index 
of refraction in dielectrics relies on mostly absorptive scattering and comes down to the original 
effective medium one [3J. 

There has been a renewed interest in systems consisting in arrays of resonators in a host medium 
since the last decade, notably owing to the emergence of metamaterials [14l|5j|. Even though many 
metamaterials now rely on non resonant unit cells ||6j I3> we abusively ignore them here when 
mentioning metamaterials. John Pendry, who first proposed metamaterial presenting electric and 
magnetic activities in electromagnetism iHIll, adopted a macroscopic point of view and devel- 
oped an approach based on field averaging to obtain the overall response of the resonators and 
the effective properties of those man made materials [31151. In this formalism the period of the 
considered medium ought to be much smaller than the free space wavelength. Nowadays those 
materials are commonly studied from the far field which precludes any knowledge of the field at 
the unit cell level, and mostly to design media with negative index [fT0] - fT2| . 

In the meantime, mainly in acoustics, there has been a large interest for materials presenting 
hybridization band gaps [[T3l - l20l . The latter, contrary to conventional photonic or phononic crys- 
tals [|2n - l24]| . present stop bands that are related to the resonant nature of the unit cell rather than 
the period of the medium. This phenomenon occurs when a resonator hybridises with the contin- 
uum of the plane waves of a homogeneous medium, giving rise to a binding and an anti-binding 
branches, that are separated by the so-called hybridization band gap [[T4ll . This effect is usually 
rather elusively justified by a level repulsion between the wave (the photon or the phonon), and 
a local resonance lfT4l . Interestingly similar structures have been proposed to decouple antennas 
||25ll26ll and are named electromagnetic band gap materials (EBG). 

At this stage, one could wonder if we are not dealing with the same physical phenomenon albeit 
named differently. Metamaterials and hybridization based materials are indeed both obtained when 
a homogeneous medium is filled with resonant elements and they both present bands of permitted 
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and prohibited propagation. There are, however, drastic differences in the common understanding 
of those phenomena. Because of the deep subwavelength nature of metamaterials, their effective 
parameters and bandwidth are usually justified by a strong near field coupling between the ele- 
ments [mm. On the contrary, hybridization effects suppose solely a coupling between free space 
waves and a single resonant element. One question arises: are those phenomena identical, and if 
so, where does dispersion come from in those materials? 

In this letter, we use a very simplified model in order to grasp the physics of an array of res- 
onators in a homogeneous medium. We first prove that, to a large extent, the response of a quasi 
ID medium is governed by a far field coupling between the individual elements. This coupling 
can be understood as Fano interferences ll27l l28l between the incoming plane waves and the field 
reemitted by the resonators. We give a phenomenological description of this effect in terms of 
the frequency response of an oscillator. We propose a simple formalism that gives the dispersion 
relation of quasi ID metamaterials using solely the far field transmission coefficient of a single 
unit cell and the period of the medium. We verify our approach and the obtained formalism on 
various media including random and wavelength scaled ones. 

Our first goal is to identify the type of coupling that can give rise to an effective negative prop- 
erty in a metamaterial. To that aim, we study numerically a very simple system consisting in an 
array of 40 split ring resonators (SRR) made of perfect electric conductor (PEC) and placed in a 
waveguide of subwavelength width. The SRR resonate around /o for transverse magnetic polar- 
ized (TM) waves and have a square shape Ao/9 wide. Such a system is known to present a negative 
permeability. We perform the simulation in a 2D domain, and hence our approach is quasi ID 
(Fig. l.a.). Because we want to separate near field contribution from far field ones, we simulate 
the structure with two types of waveguides: the first has PEC boundaries and hence no cutoff for 
the impinging TM waves, while the second has perfect magnetic conductor (PMC) boundaries 
meaning that the impinging waves are evanescent along the waveguide. This way we can select 
which coupling we allow: near field only for PMC or near and far field couplings for PEC. In 
Figure 1 .b, we map the magnetic field in the array of SRR for two different frequencies and for 
the PEC and PMC waveguides. Both show subwavelength varying patterns correspondent to a 
high effective index of refraction. More interesting are the plots of the transmission coefficients 
presented in Figure l.c for the two studied cases. For the PMC waveguide, which forbids any far 
field type of coupling between the resonators, the modes of the system are symmetrically posi- 
tioned around the resonance frequency /o of a single SRR on a narrow spectral range. This type of 
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FIG. 1. Simulated PEC (blue) or PMC (red) subwavelength waveguide filled with an array of SRR (a), 
maps of the fields for the two configurations and at two frequencies (b), transmission coefficients of the 40 
SRR for PEC (blue) and PMC (red) case (c), transmission coefficients of a single SRR for PEC (blue) and 
PMC (red) case (d). 



branch is typical of the tight-binding model commonly used in solid state physics [|29l or coupled 
resonators optical waveguides (CROWs) [|30l. and cannot explain neither the effective negative 
permeability of the array of SRR, nor any band gap. On the contrary, the transmission coefficient 
of the PEC waveguide structure displays the wanted features: it is highly asymmetric, shows a 
wide band of propagating modes below /o and a wide stop band above. The band of propagat- 
ing modes can be seen as a high permeability one, and the band gap as a negative permeability 
window. However another interpretation can be given. Above the band gap, another set of modes 
seem to be observable from the transmission coefficient. This band, well known for materials 
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presenting hybridization band gaps, is the anti-binding one, while that of the lower frequency and 
high wavenumber modes is the binding one. Those two branches are separated by a hybridization 
band gap. The response of this system can hence be explained by an hybridization between the 
resonance of a SRR and the continuum of free space waves. 

From here two consequences appear. First, in order to obtain the expected properties of an 
array of SRR, one must take into account the far field coupling between the resonators. Second, 
this system, a ID metamaterial, can be interpreted through the prism of hybridization. In order to 
obtain a deeper insight of the coupling involved, we simulate the response of a single resonator 
placed in the same waveguides. The transmission coefficients obtained are plotted in Figure l.d. 
For the PMC case, the latter displays a symmetric shape consistent with a resonant tunneling of 
the electromagnetic energy at the resonance of the SRR. When the waveguide has PEC boundary 
conditions, the transmission coefficient has a singular profile since it presents a maximum of trans- 
mission right before the resonance and a minimum slightly after: this is typical of a Fano profile 

mm. 

We prove now that Fano interferences explain the coupling between the resonators and its far 
field nature. The principle, schematised in Figure 2. a, can be explained in the following way. 
The resonator placed in the propagating waveguide scatters the incoming plane wave. Around 
resonance, the resonator accumulates energy and re-radiates waves which interfere with the un- 
scattered incoming plane waves. Now the phase of the field in the oscillator changes abruptly 
at the resonance from to tt, meaning that its contribution can interfere constructively with the 
incoming plane waves or destructively. This explains the Fano profile and gives a very interesting 
information; since the resonator experiences a phase shift at the resonance, it can "give" some ex- 
tra phase to the incoming plane wave, but it can also drastically suppress its transmission through 
destructive interferences [31]. 

The principal features of the array of SRR can be inferred quite straightforwardly from this 
approach. On one hand, before resonance, the extra phase that originates from the resonators can 
compensate for the small phase shift accumulated through propagation over a unit cell, meaning 
that the field in the array can vary on scales as small as its period, hence leading to a high effec- 
tive permeability. On the other hand, after resonance, the SRR respond in anti-phase and cancel 
the incoming plane wave which results in an hybridization band gap or equivalently an effective 
negative permeability. 

Going from one resonator to an array, the effect of Fano type interferences are cumulative, hi 
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Figure 2.b, we plot the transmission coefficients of 1, 2, 3 and 40 SRR. The number of trans- 
mission maxima and minima equals the number of resonators on each branch. Each resonator 
gives rise to a Fano interference in the spectrum created by the previous one. This has two con- 
sequences; first, it widens the band of high or negative effective permeability. Second, a system 
of A'^ resonators supports 2N modes with distinct resonant frequencies, on the binding branch 
and on the anti-binding branch. Again, this differs from a tight-binding model, which gives 
only modes symmetrically positioned around the resonant frequency /q. Naturally, an infinite 
system supports a continuum of modes resulting from a continuum of Fano interferences. Hence, 
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the coupling between the resonators can be understood simply in terms of Fano interferences of 
far field components of the spatial spectrum. 

We now utilize the formalism of the transfer matrix [|32l in order to obtain the dispersion rela- 
tion of a linear array of resonators that are coupled through propagating waves only. To do so, we 
perform calculations very similar to those done on CROWs [[30l l33l and obtain a relation of the 
formUMl: 



where k is the wavenumber in the medium, a is the period of the medium and t is the complex 
far field transmission coefficient of a single unit cell. 

This formula gives both the wavenumber k of the propagating modes supported by an array 
of resonator, but also the attenuation coefficient k in the hybridization band gaps. We first verify 
its validity by numerically studying various types of resonators. Our simulations are lossless and 
all dimensions are given relatively to the wavelength at the first resonance of the considered unit 
cell Xq. We study the response to a plane wave of an array of 40 resonators in a waveguide. We 
simulate first a linear array of subwavelength SRR (Fig. [3ja). Then we apply our approach to 
Mie resonators, that is, squares of n=10 dielectric material (Fig. 3.b). We also simulate using 
acoustic waves the soda cans, that is the Helmholtz resonators, used in ll35ll . this time using a 3D 
waveguide (Fig. 3.c). Finally we study Ao/4 high, infinitely long slits in PEC, a building block of 
designer's plasmons [l36l[37ll (Fig. 3.d). We calculate the spatial spectrum of the electromagnetic 
or acoustic fields in those structures at each frequency and map the dispersion relations Figure 3.a- 
d. Meanwhile, we simulate each resonator independently in a waveguide and obtain its complex 
far field transmission coefficient. We superpose on the maps our theoretical dispersion relation 
and plot the calculated attenuation coefficients in the gaps using equation [TJ 

Figure |3] shows that for every medium under study, our simple formalism, based on propagating 
waves only, exactly reproduces the dispersion relation obtained by rigorously simulating the entire 
structure. This means that those structures, although organized on a deep subwavelength scale, are 
coupled mostly through propagating waves. Furthermore, one each structure gives rise to a wide 
and efficient hybridization band gap. For each medium, the frequency /o separates the propagating 
modes and the hybridization band gap, which justifies the phenomenological explanation given 
above OTl . This means that all those quasi ID materials have in common a unique physical basis: 
an hybridization between the resonator and the propagating waves. This approach permits to unify 
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FIG. 3. Structures and parameters, maps of the numerically obtained dispersion relations with superposed 
our theoretical predictions, and calculated attenuation coefficients for the: TM excited Aq/S spaced array 
of SRR (a), TM excited Aq/V spaced array of Mie resonators (b), acoustically excited Ao/12 spaced array 
of soda can or Helmholtz resonators (c), TM excited Ao/20 spaced array of PEC slits, that is, a designer's 
plasmons surface (d). 

several distinct ones. For instance designer's plasmons [|36]438l, which are all based on a resonant 
unit cell, can be analysed under the same prism than negative permeability [HI, permittivity [[39] - 
[4T| or Young's modulus materials [[T3l[42ll43l . This is in agreement with the fact that an interface 
between a negative and a positive effective property medium supports surface waves [[40l[44]| which 
we utilized in order to beat the diffraction limit from the far field [l35ll45] - [47l . 

From these results, it is clear that near field couplings play a very little role in the dispersion 
observed in the systems studied here. This leads us to ask the following question: does one need to 
fulfill the effective medium condition of a deep subwavelength period in order to design media with 
interesting properties? To answer this question and demonstrate the possibilities offered by this 
Fano interference type of coupling, we present additional simulations. We want to demonstrate that 
hybridization band gaps, or equivalently bands of negative effective properties, exist even when 
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FIG. 4. Structures and parameters, maps of the numerically obtained dispersion relations with superposed 
our theoretical predictions, and calculated attenuation coefficients for the: TM excited Ao/2.4 spaced array 
of SRR (a), TM excited Xq/U spaced array of Mie resonators (b), TM excited and randomly positioned 
array of SRR (c), TM excited Aq spaced array of Mie resonators (d). 

the spacing between the resonators is increased and resist to disorder. To that aim we simulate four 
more systems. The first one consists in an array of SRR separated by Ao/2.4, a distance comparable 
to the period of typical Bragg crystals (Fig. |4}a). Similarly we study a Aq/I.V periodic array of 
dielectric squares of index n=10 (Fig. 4.b). We simulate a chain of randomly positioned SRR 
(average period a of Aq/S, Fig. 4.c). Finally we study an array of Mie resonators of period of 
Ao (Fig. 4.d). Again we use equation [T] to superpose our theoretical dispersion relations and plot 
accordingly the attenuation coefficient in the band gaps. Those curves totally fit the simulations. 

Figure 4 shows a very important feature common to all studied media: there exist in all the 
dispersion relations bands of propagation and hybridization band gaps regardless of the medium's 
order or period, thanks to the far field coupling between the resonators. Their positions, however, 
can be shifted due to the combined effects of interferences in between the unit cells and the reso- 
nant nature of the latter. Nevertheless, this proves that material made of resonant unit cells, even 
spaced by distances of the order or larger than the wavelength or random, can present negative 
effective properties. 
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Figure 4 also proves that Bragg bands and band gaps coexist with hybridization dilute systems, 
as was already proved in IHHIHl, hence providing an even greater flexibility for designing ma- 
terials with unusual properties. Due to their resonant nature, hybridization band gaps are much 
more efficient than Bragg ones. Mixing Bragg and hybridization effects could possibly lead to 
very interesting new phenomena. 

We note that disorder plays a different role for Bragg and hybridization types of bands. Indeed, 
in Figure 4.c, which shows the results of an array of SRR randomly positioned, the Bragg band gap 
has disappeared, while the hybridization one remains, a fact demonstrated in the past in acoustics 
[fT4l - [T8l . This underhnes the potential applications of hybridization band gaps in the visible. 

This work finally appeals some concluding remarks. We have shown here that dispersion in 
arrays of resonators can be understood in terms of a far field type of coupling between the elements 
that arise from Fano interferences between the incoming waves and the waves re -radiated by the 
resonators. This approach gives a physical mechanism to the hybridization effect often referred 
to in acoustics and applies to metamaterials, locally resonant materials, designers plasmons and 
EBG structures. It demonstrates that wide bandwidth negative effective properties or hybridization 
gaps do not require any near field coupling and exist in random and wavelength scaled media. We 
believe that mixing Bragg and hybridization effects makes the physics of man made materials 
even richer and their properties more flexible. We underline however that some materials based 
on resonant unit cells can be also coupled through near field interactions, such as inductive or 
capacitive couplings in electromagnetics. This is typically the case for closely spaced SRR in 2D 
or 3D which can support magneto-inductive waves [48 J, or bubbles in acoustics [TTSl |49l . Such 
near field coupling may be taken into account in our formalism, but this is out of the scope of 
this paper. Our approach requires a generalization to 2D and 3D materials. This could permit to 
design materials with exotic properties quite simply. Initial results on our experimental system of 
ll35ll support our intuition that 2D or 3D systems behave exactly as the ID one studied here. This 
suggests that since the Fano based dispersion in those media is linked to the forward scattering 
cross-section of the individual resonators, metamaterials are inherently non-local. Finally this 
idea should be transposable to all systems that support waves and contain resonant inclusions, 
which are very common in nature. 

The authors acknowledge fundings from Orange Labs, and F. Lemoult acknowledges funding 
from French "Direction Generale de I'Armement". 
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